MeV-GeV emission from neutron-loaded short gamma-ray burst 

jets 



Soebur Razzaque and Peter Meszaros 

Department of Astronomy & Astrophysics, Department of Physics, Pennsylvania State 

University, University Park, PA 16802 
soeb@astro.psu.edu, nnp@astro.psu.edu 

ABSTRACT 

Recent discovery of the afterglow emission from short gamma-ray bursts sug- 
gests that binary neutron star or black hole-neutron star binary mergers are the 
likely progenitors of these short bursts. The accretion of neutron star material 
and its subsequent ejection by the central engine implies a neutron-rich outflow. 
We consider here a neutron-rich relativistic jet model of short bursts, and in- 
vestigate the high energy neutrino and photon emission as neutrons and protons 
decouple from each other. We find that upcoming neutrino telescopes are unlik- 
ley to detect the 50 GeV neutrinos expected in this model. For bursts at z ~ 0.1, 
we find that GLAST and ground-based Cherenkov telescopes should be able to 
detect prompt 100 MeV and 100 GeV photon signatures, respectively, which may 
help test the neutron star merger progenitor identification. 

Subject headings: gamma rays: bursts — gamma rays: theory — ISM: jets and 
outflows — radiation mechanisms: non thermal — stars: neutron 

1. Introduction 

The nature and origin of short gamma-ray bursts (SGRBs), with t 7 < 2 s duration, 
were largely conjectural until recently. This changed dramatically with the detection of the 
afterglow and the identification of the host galaxy of GRB 050509b at a redshift z = 0.226 
(Gehrels, et al. 2005), as well as subsequent observations of other short burst afterglows, 
e.g., GRB 050709 (Fox, et al. 2005; Villasenor, et al. 2005), GRB 050724 (Barthelmy et al. 
2005), etc. These observations revealed that the short bursts are located mainly in elliptical 
hosts, while some are in spirals or irregulars, as expected from an old population such as 
NS-NS or NS-BH mergers. They also provided, for the first time, details about the X-ray, 
and in some cases optical and radio afterglow lightcurves, including in two cases evidence 
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for a jet break, qualitatively similar to the relativistic jet afterglows of long GRBs [see, e.g., 
Zhang & Meszaros (2004) for reviews of long GRBs]. The average redshifts are lower than for 
long bursts, indicating an isotropic-equivalent energy of SGRBs approximately two orders of 
magnitude lower than for long ones. 

We investigate here a model of SGRB jet which is initially loaded with neutrons and 
fewer protons, as expected from a neutron star composition. High energy emission in this 
model is then expected from the photosphere of the relativistic jet outflow, including both 
thermal and non-thermal components, in addition to the radiation expected outside of it. 
We describe our jet model in Section 2, high energy emission processes in Section 3, and 
expected neutrino and photon fluxes in Section 4 and 5 respectively. We summarize and 
discuss our results in Section 6. 



2. Neutron loaded jet dynamics 

We assume the total isotropic-equivalent energy outflow in the SGRB jet of L = 10 50 L5o 
erg/s. The jet is loaded with protons and neutrons with an inital neutron to proton number 
density ratio £ G = n' n /n' p = 10£ O) i in the comoving frame. 1 The total mass outflow rate in 
such a neutron-rich jet is M = 4irr 2 c(l + ^ )n' p m p initially, since m p ~ m n m e and the 
thermal energy is negligible. A useful quantity is the usual total energy to mass flow ratio 
7] = L/Mc 2 , characterizing the dimensionless entropy or the baryon loading of the jet. In 
a typical long GRB (t 1 > 2 s) the final bulk Lorentz factor of the 7-ray emitting fireball is 
equivalent to rj. In the case of SGRBs, we assume 77 = 316t7 2 . 5 , similar to that of the long 
GRBs, motivated by observation. 

The outflow starts at a radius R Q = 10 6 i? Oi6 cm which is a few times the Schwarzschild 
radius r g = IGM^fc 2 of a solar mass black hole created by binary mergers. The inital 
temperature of the plasma outflow is 
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where a = 7r 2 /c 4 / '15 (he) 3 = 7.6 x 10 -15 erg cm -3 K~ 4 is the radiation density constant. The 
optically thick hot plasma expands adiabatically due to the radiation pressure. The comoving 
temperature drops as T'(r) oc R /r and the plasma expands with an increasing bulk Lorentz 
factor r(r) oc r/R a with the radius r following the adiabatic law. Electrons (both e ± and 



1 We use primed variables in the relativistic plasma frame and un-primed variables in the observer frame 
or laboratory frame 
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n e ~ n p ) are coupled to the photons via Compton scattering, and protons, which are coupled 
to the electrons via Coulomb interaction, are held together in the expanding optically thick 
plasma. Neutrons, however, are somewhat loosely coupled to the protons by elastic nuclear 
scattering with a cross-section a np »3x 10~ 26 cm 2 , which is roughly a fraction 1/20 of the 
Thomson cross-section (orh ~ 6.65 x 10~ 25 cm 2 ). This changes considerably the dynamics of 
a neutron- rich jet [e.g. Derishev, Kocharovsky & Kocharovsky (1999); Bahcall & Meszaros 
(2000); Beloborodov (2003)] as we discuss next. 



r 



The protons and neutrons are coupled together with a common bulk Lorentz factor 
r) ~ Tp(^) — r(r) as long as the n-p collision time t' np ~ {n' p a np c)~ l is shorter than 
the plasma expansion time t' exp ~ r/cY{r). The critical dimensionless entropy for which n-p 
decoupling happens, from the condition t' n , 
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As a result, the neutrons and protons in our SGRB jet model (n > r] np ) decouple at a radius 
where the nuclear scattering optical depth r' np ~ n' p a np R np /T(R np ) = 1 as 
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The bulk Lorentz factor of the neutrons does not increase further and saturates to a final 
value of 
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at r = R np . The neutrons start to lag behind the protons at r > R np , and the total energy 
outflow due to protons is 



L — L — Y n jMc z i = L 1 
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7 x 10 49 ergs/s, 



(5) 



which is close to the initial value of L since T n j < n. Here M ~ M/(l + is the 
proton mass outflow rate after decoupling. The outflow is still optically thick at r = R 



up 



since a^h — 20a np . The decoupling of radiation from plasma electrons takes place for a 
critical value for the dimensionless entropy, for which the Thomson scattering time scale 
t' Th ~ (n'gOThc) -1 is equal to the plasma expansion time, is given by 
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for our SGRB parameter L 49 , 8 — L/(7 x 10 49 ergs/s). Note that the dimensionless entropy 
of the proton outflow fj = L/Mc 2 is related to 77 ra d by 
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In this case (fj > ?) ra d) the outflow becomes transparent during the acceleration phase T oc r, 
at a radius r = R T where the Thomson optical depth r^ h ~ n' e R T /T(R T ) = 1, 
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which is close to i? np . The final Lorentz factor of the (proton) outflow is given by (Rossi, 
Beloborodov & Rees 2005) 
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which is achieved at a radius r ~ R T for this 77. This suggests a possible way to generate 
the high bulk Lorentz factor estimated in short GRBs without any spectral lags in the data 
(Norris & Bonnell 2006). 

Note that for 77 < i] np , one has Y n j = F p j = r\ for a baryon loaded jet outflow satisfying 
the condition r\ < i] ra d- Here ?7 ra d is the initial critical radiation entropy for combined proton 
and neutron outflow, defined with an equation similar to equation (6) as 
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We have plotted the final proton and neutron bulk Lorentz factors as functions of the 
baryon-loading parameter in Figure 1. For a given total energy outflow rate L, the maximum 
allowed value of r], from equation (10), are shown as three vertical solid lines for three values 
of initial neutron to proton ratio £ G =1, 5 and 10. Note that for r\ ~ 3I6772.5, a typical 
value used for long GRB modelling, the neutron and proton components always decouple 
[see equation (2)] in case of SGRBs because of a lower L. On the contrary, n-p decouple 
only for £ G > 4 in case of a long GRB with rj ~ 3I6772.5, L = 10 52 L 52 erg/s and R a = 10 7 R o j 
cm. Also note the high T p j values of 519, 572 and 624 for £ = 1, 5 and 10 respectively in 
Figure 1 from equation (9). 



3. Inelastic neutron-proton scattering 



The elastic component of the total np cross-section, which dominates at lower energy, 
drops rapidly above the pion production threshold of as 140 MeV, where the inelastic np 
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scattering cross-section is a n ~ a np . Hence, the condition n > i] np (and subsequently fj > 
^rad) for n-p decoupling implies that there will be pion production by the threshold processes, 
with roughly equal (total unit) probability to produce a 7r ± pair or (and) 2n°, 



nnn + ; n + — > fj+v^ — > e + u e u ll u IJl 



pn -> <] ppTT ; tt v^->e v e v^v^ (11) 

pn7r° ; 7T° — > 77 

The secondary pions produced by the inelastic np scattering are at rest in the center-of-mass 
(cm.) frame of interaction. However, the neutrons are not cold in the comoving outflow 
frame as we discuss next, and the cm. frame acquires a velocity relative to this frame. 

The Lorentz factor of the neutrons relative to the outflowing protons, at the radius 
r ~ R T ~ R np , is given by 



r rel = r p r n (i - (3 p( 3 n ) ~ U^L + ^) « 2.8 
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Thus, in (proton) outflow rest frame the neutron energy is e' n = Y TC \m n c 2 and the Lorentz 
factor of the cm. is given by 7^ m = (e' n + m p c 2 ) j '(m^c 4 + m^c 4 + 2e' n m p c 2 ) 1 ^ 2 xs 1.4. The 
energies of the decay neutrinos are ~ ~ 307^ m MeV from -k ± , e' Ue ~ e' Pe ~ 307^ m 
MeV and e' v ~ ~ 507^. m MeV from /i 1 * 1 . The 7r° decay photon energy is ~ 707^. m 
MeV. The observed energies of the decay products would be 

{26 GeV ; [y^ and u e from ^ and ji^ 1 ) 

43 GeV ; (^ from //±) x r P)2 . 8 . (13) 

60 GeV ; (7) 

Here r p>2 . 8 = r P)/ /624. 

Note that additional pions may be produced via pp — > pn7r + and nn — > np7r~ inelastic 
interactions with similar cross-sections as the inelastic np interaction. The nn process in 
particular may be important for neutron destruction in the outflow (Rossi, Beloborodov & 
Rees 2005). Neutrinos and photons produced from nn and pp processes, however, are of 
much lower energy as the relative velocities between the n-n or p-p components (e.g., arising 
from a variable outflow) may not exceed that between n-p in equation (12). 



4. Neutrino emission 

The inelastic n-p scattering opacity is ~ r' np ~ 1 at r ~ R np ~ R T and each proton 
produces, on average, a 7r ± pair half of the time. The decay neutrinos escape the plasma 
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outflow as soon as they are created. The neutrino and anti-neutrino flux rate from an 
SGRB at z ~ 0.1 (corresponding to a luminosity distance Dl ~ 10 27 D 2 7 cm for a standard 
cosmology with Hubble constant H ~ 70 km s _1 Mpc -1 ) is 

K - K - - ™* * 47tD J - , « 5.7 x 10- 6 -^^ cm- 2 s" 1 . (14) 

4:71 JJ L L pjlTlpC D 27 L Pj2 .8 

The neutrino-nucleon (iMT) total cross-section for detection is cr„jv ~ 6 x 10~ 37 cm 2 at 
ej, 60 GeV [see, e.g., Table 1. in Gandhi, et al. (1998)]. Upcoming kilometer scale water 
Cherenkov detectors such as ANTARES and KMSNeT are designed to have sensitivity for 
detecting neutrinos of energy > 50 GeV. The total number of nucleons in a gigaton detector 
is Njy ~ 1.2 x 10 39 , which folded with the cross-section and muon neutrino flux from a typical 
SGRB at z ~ 0.1 yields an event rate of 

Nfa + event) ~ 2N^a uN N N « 0.01 ^g.s^^i.s s^Gton" 1 (15) 

per burst. Hence it is unlikely that upcoming neutrino telescopes would be able to detect 
these neutrinos, based on our model as dicussed here, from individual SGRBs. The number 
of events from diffuse flux can be calculated assuming a burst rate of ~ 300/47T yr _1 sr _1 
(approximately 30% of the long GRB rate) distributed over the whole sky and an optimistic 
27T sr angular sensitivity of the detector as (300/47r) x 2ir x iV(z/ M + event)t « 1 yr _1 . 
Here, we have assumed t ~ 1 s as the typical duration of the SGRBs. 



5. Photon emission 

The 7-ray emission from an SGRB, however, is more promising. Although the plasma 
outflow becomes transparent to Compton scattering at a radius r ~ R T ~ R np ~ 10 s cm, 
not all the high energy photons due to 7r° decay can escape. Only those created below a skin 
depth of R T escape with a probability (Bahcall & Meszaros 2000) 

P n o ~ r' np {R T ) ~ R np /R T « 0.4. (16) 

The corresponding number flux of 60 GeV 7r°-decay photons at Earth is 



N, n o ~ - J^ L « 2.0 x 10- 6 - j I 9 ' 8 cm" 2 s" 1 , (17) 

47T J D 2 r pJ m p C 2 ^27^,2.8 



which is below the detection threshold of GLAST "Large Area Telescope" (LAT), whose 
effective area is ~ 10 4 cm 2 in this energy range (Gehrels & Michelson 1999). However, 
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ground-based Cherenkov telescopes with < 100 GeV threshold may detect these photons. 
In particular, Milagro has an effective area of ~ 5 x 10 5 cm 2 for a zenith angle < 15° in this 
energy range (Dingus 2004). For a yearly SGRB rate of ~ 300/yr, this implies a possible 
detection rate of < 5 yr _1 , assuming a 90% duty cycle. 

The high energy photons from 7r°-decay created below R T interact with thermal photons 
of energy e' Jt ~ T' (R /R T ) 3.6 keV, which are carried along with the plasma. The 
number density of these photons is high, and all 7r°-decay photons of energy e 7i 100 
MeV originating below the r^ h ~ 1 skin depth produce pairs satisfying the condition 
e 7 t e '~ f i ~ rn 2 c A . Assuming equal energy distribution among final particles, charged pion decay 
and 77 — > e ± have a characteristic energy ~ 357^ m MeV, and initiate electromagnetic 
cascades by inverse-Compton (IC) scattering of thermal photons, creating new generations 
of pairs which interact again with the thermal photons. Contrary to the optically thick 
case of r^ h 3> 1, where all upscattered photons are absorbed by thermal photons producing 
e ± pairs, these cascades are not saturated (Derishev, Kocharovsky & Kocharovsky 1999). 
The upscattered photons eventually escape from the photosphere with characteristic energy 
e' c ~ m e c 2 below the e ± pair production threshold energy. 

The fraction of the 7r ± and 7r° energy carried by the escaping 7-rays depends on the 
details. An order of magnitude estimate for this fraction is £ 7 = 0.1£ 7i _i (Rossi, Beloborodov 
& Rees 2005). Thus, the number flux of the escaping 7-rays of energy 

e 7 , c ~ T pJ m e c 2 « 320 r p , 2 . 8 MeV (18) 

at Earth is 

Kc - ttJ^ 2 ( ^-^^-) - 8-6 x 10- 4 ^#^ cm" 2 s" 1 . (19) 

The escaping Comptonized photon spectrum of the cascade is d/V 7 /de 7 oc e~ q , where q « 5/3 
for e 7 < e 7jC and q ~ 2.2 for e 7)C < e 7 < e 7)inax (Belyanin and Derishev 2001). The 
maximum photon energy is e 7jmax ~ l'e e j,t^pj ~ 20 GeV. The fluence threshold for GLAST 
is ~ 4 x 10 -8 ergs cm -2 for a short integration time (Gehrels & Michelson 1999). Thus the 
Comptonized photon fluence e^ >c F €j c ~ e 7iC iV 7iC t w 4.4 x 10~ 7 ergs cm -2 for a typical SGRB 
of t ~ 1 s duration at z ~ 0.1 should be detectable by GLAST (see Figure 2). 

Apart from the 7r°-decay photons of energy 60 GeV and e ± pair cascade photons of 
energy 320 MeV, the epectrum of radiation from the photosphere of an SGRB would also 
contain thermal photons of peak energy e 7i & ph 2.82e 7 t r p j 6.4r Pi2 . 8 MeV (see Figure 2). 
The corresponding number flux on Earth of these photospheric photons of peak energy is 

= i^r " a5 cm " 2 (20) 

47T7J L e 7)fe D 27 L P)2 .8 
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The photon fluxes at various energies discussed above are associated with the jet pho- 
tosphere, and last as long as the usual burst prompt 7-ray emission. Hence, these emissions 
would be contemporaneous with the emission at ~ MeV associated with internal shocks or 
other dissipation prcesses above the photopshere, aside from a small time lag between the 
onset of the two events of the order of a millisecond. We discuss emission from internal 
shocks next. 



The internal shock model of plasma shells colliding at some radius above the photosphere 
is the most commonly discussed mechanism to produce the observed 7-rays in the fireball 
shock model. The average bulk Lorentz factor of the colliding shells is Tj ~ Y p j 6241^2.8- 
With a millisecond variability time t v = 10 _3 t„ 5 -_3 s, the internal shocks take place at a 
radius n ~ 2T}ct v « 2.3 x 10 13 cm. Note that this is much smaller than the /3-decay radius 
~ 3.4 x 10 15 cm. Hence, most neutrons have not had time to decay until well beyond the 
internal shock radius. 

The equipartition magnetic field in the shock region is B = (2eBL/rfc) l l 2 ~ 9.1 x 10 5 G 
with an equipartition fraction Sb = O.le^-i. The shock-accelerated electrons are assumed 
to follow a power-law energy distribution dN' e /d-y' e oc r )'~ a for 7^ > n f' e min — e e m p /m e with 
a ~ 2.4, which results in synchrotron radiation with a peak energy 



Here e e = 0.1e e -i is the fraction of fireball's kinetic energy transferred to the electrons by 
shocks. The corresponding number flux at Earth of these peak synchrotron photons is 



The typical synchrotron radiation spectrum from power-law distributed electron energy is 
oc €j ~ e 2 2 for €j > c-y^m (sgg Figure 2). The observed typical GRB spectra 

dN 1 /de 1 oc e^ 1 , is slightly different from the theoretical synchrotron spectrum with an index 
-3/2. 

We have plotted the energy spectra of synchroton radiation from the internal shocks, 
and thermal and IC cascades from the jet photosphere in Figure 2. We used L = 10 50 L 50 
erg/s, R Q = 10 6 cm and 77 = 316 for all curves and £ = 1, 5 and 10 for the dashed, dotted 



5.1. Comparison with internal shock emission 
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and dotted-dash curves respectively. The peak energy of thermal photons (e 7j &) decreases 
from 6.4 keV for £ G = 10 to 3 keV for £ G = 1 as the photosphere radius (R T ) in equation 
(8) increases with decreasing £ G . The variation of the peak energy (e 7jC ) of the photospheric 
cascade is directly proportional to the variation of F p j with £ Q and increases with £ G . We 
have normalized these spectra by integrating the corresponding differential number fluxes 
from e 7jC up to the respective maximum energy (e 7jmax ) and equating to the respective total 
number flux in equation (19) with fixed £ 7 = 0.1, the fraction of the total energy carried 
by the pions. The variation of the synchrotron peak energy (e 7jm ), in the range 330-400 
keV, with £ D can be explained similarly. We have normalized the synchrotron spectra by 
integrating the corresponding differential number fluxes from e 7jTn up to the energy 10e 7)m ~ 
few MeV, the typical energy range for the BATSE and BAT detectors, and equating to 
the respective total number flux in equation (22). We have plotted an extension of the 
synchrotron spectra up to 1 GeV, in order to compare them with other spectral components. 

The synchrotron spectrum does not extend very far above e 7 , m , for the parameters used 
here, from a calculation of the maximum electron Lorentz factor. The energy at which 
the Rayleigh- Jeans portion of the thermal photospheric photons start to dominate over the 
synchrotron spectrum is e 7 ~ 2 MeV for £ G = 5 and 10 curves in Figure 1. (For £ = 1, 
the thermal component is small over the whole energy range.) Above this, besides the 
Comptonized photospheric cascade spectrum, there could also be an IC scattered component 
of the synchrotron photons (SSC), whose importance is characterized by the parameter 
Y = [—1 + (1 + Ae e /e B ) l / 2 }/2 0.6 for our present model. The corresponding minimum IC 
peak is 7^ min ~ 4 x 10 4 higher than the synchrotron peak, at a photon energy e 7jmj ic ~ 10 
GeV. Hence an SSC component is not important in this case (Y < 1). Thus, at energies above 
which the sycnhrotron spectrum cuts off, the Comptonized photospheric cascade photon 
spectrum would be clearly distinguishable by GLAST (see Figure 1). Although for some 
choice of parameters an SSC contribution could become significant at energies > 10 GeV, 
this is expected to be at energies much higher than the peak of the Comptonized cascade peak 
e 7jC ~ 300 MeV. The maximum SSC photon energy, in the Thomson limit 7g <C m e c 2 /e 7)m , 
would be e 7j ic <sC ^1^1/ e -y,m ~ 300 GeV. At energies well above ~ 10 GeV the two 
components (Comptonized photospheric cascade and SSC) may be hard to distinguish, but 
in the range of ~5 MeV - 10 GeV the Comptonized photospheric cascade should dominate. 
Below the peak e 7 < e 7;C , the Comptonized photospheric cascade component is lower than the 
thermal photospheric peak, and at even lower energies it is also lower than the synchrotron 
component. 



6. Discussion 



If SGRBs are produced as a result of double neutron star or neutron star-black hole 
binary mergers, the neutron-rich outflow may result naturally in a high proton Lorentz factor 
T p j > 600. Besides the usual non-thermal (e.g. internal shock) synchrotron spectrum, this 
implies the presence of a thermal photospheric component peaking around ~ 10 MeV. The 
n-p decoupling will lead to higher energy photons (> 0.1 GeV) from neutral and charged 
pion decay cascades, and neutrinos (> 30 GeV) directly from pion and muon decays. There 
may be additional X-ray emission when the neutron-decay proton shell with initial lower 
Lorentz factor impact with the decelerating proton shell initially moving with a high Lorentz 
factor (Dermer & Atoyan 2006). The afterglow light curves can also be affected by the 
neutron decay (Derishev, Kocharovsky & Kocharovsky 1999; Rossi, Beloborodov & Rees 
2005). Here, we have concentrated on the prompt (t < 2 s) high energy (MeV to tens 
of GeV) signatures of neutron-rich outflows. For short bursts at z ~ 0.1 the ~ 50 GeV 
neutrinos are unlikely to be detectable with currently planned Gigaton Cherenkov detectors. 
However, the pion decay photons at ~ 60 GeV escaping from whithin a r^ h ~ 1 skin-depth 
of the photosphere should be detectable with high area, high duty cycle Cherenkov detectors 
such as Milagro, at an estimated rate of 5 bursts per year. The pion decay photons created 
below the photosphere will result in a Comptonized cascade spectrum which, peaking in the 
range > 0.3 GeV, dominates the non-thermal shock synchrotron-IC spectrum, and should 
be detectable with GLAST. Both the prompt GeV photon components discussed here are 
essentially contemporaneous with the usual MeV range short GRB prompt emission, t 7 < 2 
s. Their detection or lack thereof could constrain the neutron content of the outflow and the 
nature of the short burst progenitors. 

We thank the referee for comments which have helped us to improve the paper. Research 
supported in part through NSF AST 0307376 and NASA NAG5-13286. 
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Fig. 1. — The final bulk Lorentz factor of the proton (r p /) and neutron (T n ,f) components, 
using equations (4) and (9) respectively, in the GRB jet outflow as functions of the baryon- 
loading parameter (77). The dashed, dotted and dotted-dash curves corresponds to the initial 
neutron to proton ratio £ c = n' n /n' — 1, 5 and 10 respectively. The other parameters used 
are L = 10 50 L 50 erg/s and R Q = 10 6 cm for all curves. For rj < r] np , defined in equation (2), 
T p j = Y n j = i] as shown by the single line. The protons and neutrons are held together in 
the outflow with leptons and photons in this case. For r\ > i] np the protons and neutrons 
decouple as shown by the splitting of the curves at r\ = i] np = 148, 172 and 227, and 
T p j > r n j. Note that the n-p decoupling takes place for lower rj as £ G increases with other 
parameters fixed. The vertical lines at rj = 322, 374 and 493 correspond to the limiting value 
of 77 = ?7 rad defined in equation (10). 
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Fig. 2. — Photon spectra of the expected photospheric signals (thermal and IC cascade 
by pion decay electrons and photons) and the synchrotron spectrum from internal shocks, 
compared to the GLAST LAT sensitivity. The dashed, dotted and dotted-dash curves cor- 
responds to the initial neutron to proton ratio £ Q = n' n /n' p — 1, 5 and 10 respectively as in 
Figure 1. The other parameters used are L = 10 50 L 50 erg/s, R Q = 10 6 cm and rj = 316 for 
all curves. To calculate the internal shock synchrotron spectra we used a fixed variability 
time t v = 10~ 3 s. The peaks of the spectra are correlated with the final bulk Lorentz factor 



of the proton outflow: Y p j 
text for more details). 



519, 572 and 624 respectively for £ Q = 1,5 and 10 (see main 



